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Taking into account both the scattering and the velocity saturation of carriers, we examine the 
" shallow- water" instability of the two-dimensional electron gas in a field effect transistor. It is shown 
that both the scattering (which is analogous to friction in a shallow-water channel) and the carrier 
velocity saturation lead to damping of the plasma wave instability. Threshold diagram of instability 
is calculated. The actual device parameters required for observation of plasma wave generation are 
compared with those reported in recent sub-terahertz emission experiments. 
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Recently, the interest in the plasma wave generation 
mechanism [1] in an ultrashort-channel Field Effect Tran- 
sistor (FET) has resumed. As it was demonstrated 
in Rcf.[l], the behavior of high-density two-dimensional 
electron gas (2DEG) can be described by equations anal- 
ogous to the hydrodynamic equations for shallow water, 
the plasma waves being similar to waves in this medium. 
With the carrier scattering neglected, it is shown that 
asymmetric boundary conditions with a fixed voltage at 
the source and fixed current at the drain lead to an insta- 
bility of steady state with a dc current. The FET chan- 
nel can be regarded as a resonator, with a plasma- wave 
generation mechanism similar to that known for self- 
excitation of jets and organ pipes. Both the cigenfrequcn- 
cies and the instability increment have been found for 
zero-dissipation electrons. In Refs.[2,3] a strongly non- 
linear instability mode, in which shock waves are formed 
in the channel, was analyzed. In Ref. [4], the instability 
increment and threshold were calculated for the realistic 
case of nonzero carrier scattering at arbitrary currents 
and electron momentum relaxation times. The scatter- 
ing results in the narrowing of the instability region, with 
the increment decreasing. The instability vanishes at a 
certain critical magnitude of friction. Computer simu- 
lations [4] provide a strong support for theoretical pre- 
dictions. The plasma wave instability in short-channel 
FETs is very important for application of these devices as 
high-power sources and detectors in the terahertz (THz) 
frequency range[5]-[10](see also [11]). The sub-THz emis- 
sion recently observed in short-channel FETs [12,13] was 
attributed to the plasma wave generation mechanism [1] . 

In the present paper, we analyze the plasma wave in- 
stability in the presence of both the scattering and the ve- 
locity saturation of carriers. We plot the instability dia- 
gram for device parameters required for plasma wave gen- 
eration. We demonstrate that the 2DEG samples studied 
in [12,13] have an insufficient carrier mobility and, hence, 
the observed sub-THz emission cannot be understood in 
terms of the plasma wave generation scenario. 



Steady state 

According to Ref. [1], the high-density 2DEG in an 
ultra-short channel FET can be described by the follow- 
ing equations analogous to the hydrodynamic equations 
for shallow water: 
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Here, the voltage swing, U = U gc — Ut, corresponds to a 
shallow water level, U gc is the local gate-to-channel volt- 
age, Ut the threshold voltage, V the local electron flux 
velocity, m the effective mass, and r the momentum re- 
laxation time due to collisions of electrons with phonons 
and (or) impurities. Eq.(l) is the continuity equation, 
in which the Shockley relation for the gradual channel 
approximation, ./V = OU/e, is taken into account. Here, 
TV is the surface electron concentration, and C is the ca- 
pacitance per unit area. Note that there is no direct hy- 
drodynamic analogy for the friction term ^ in the Euler 
equation specified by Eq.(2). 

It has been shown [1] that in the absence of 
friction(r — > oo), which is further named the "clean 
limit," the steady electron flow with a constant U, V is 
unstable against small perturbations 6V, SU ~ exp(— iut) 
under the following boundary conditions: 



u(o,t) = u s 



U(l,t)V(l,t)=j/WC, 



(3) 



where U s is the fixed source(a; = 0) potential, j the cur- 
rent fixed at the drain (x = I), and W the gate width. 
In the steady state, U = U s , V = V s = j/(WCU s ). The 
real and imaginary parts of u>{ see Figs.l, 2) are given by 
[1] 
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where M — V/S plays the role of the Mach number in 
hydrodynamics, S = y/ eU/m is the local plasma wave 
velocity, and n is an odd integer for \M\ < 1 and even 
integer for \M\ > 1. Actually, the instability is related 
to the plasma wave amplitude enhancement caused by 
certain boundary reflections [1]. Two different plasma 
waves with wave vectors = ±w/ (S(1±M)) propagate 
down- and upstream, respectively. The argument of the 
logarithm in Eq.(4) contains the round-trip gain factor, 
(1 + M)/(l — ill), which is the product of the reflection 
coefficients at two boundaries. The plasma wave propa- 
gates downstream during the time, l/(S(l+M)), whereas 
the back propagation takes a longer time, l/(S(l — M)). 
When M — > 1, the round-trip time of the plasma wave be- 
comes infinite, and, therefore, both the proper frequency 
and the instability increment vanish ( see Figs. 1,2 ). 

The friction of electron flow results in a spatial de- 
pendence [14] of both the voltage and the velocity. In 
this case, the instability problem for small perturbations 
superimposed on the steady state flow V, U becomes ex- 
tremely difficult. Fortunately, this problem can be solved 
by the high-order mode method suggested in Rcf.[4]. 
This approach makes it possible to deduce both the insta- 
bility increment and the threshold for arbitrary currents 
and friction strengths. As expected, the scattering leads 
to instability damping. At a certain critical friction, the 
instability vanishes. 

The special interest of the present paper is in the fric- 
tion term which has always been assumed to contain a 
constant momentum relaxation time r previously [14,4]. 
We argue that the drift velocity in real 2D systems at 
high electric fields approaches a certain saturation value 
V sa t- In the conventional FETs just the drift velocity 
saturation (DVS) mechanism determines the saturation 
current. In order to account for the possible influence of 
the DVS effect on the instability, we further investigate 
the realistic model with the momentum relaxation time 
in Eq.(2) replaced as 



(5) 



1 + iiE/V sat ' 



where fi is the low-field mobility of the 2D system, and 
E = — 4^- is the electric field. We further demonstrate 

ax 

that the DVS effect leads to an extra damping of the 
plasma wave instability. 

Using Eqs.(l,2), we obtain for the steady state 



source. Then, we introduce the dimcnsionlcss parame- 
ter p = associated with the DVS effect. Previous 
studies [14,4] correspond to the p = case. 

We are further interested only in the case of a subsonic 
flow case, when < M < 1. According to Eq.(6), at finite 
7 the Mach number increases from the initial value M s at 
the source to some drain value Md = M(l). Note that at 
fixed p the increasing solution to Eq.(6) is possible when 
(f>(M„p) > or M s < M™ ax {p) = N /p 2 /4+l-p/2. If 
M s < M™ ax , the electron flow velocity saturates at the 
drain, where the electric field is the strongest (i.e., E ~ 
— > oo ), as the expression in the denominator in the 
right-hand part of Eq.(6) goes to zero. The Mach number 
at the drain cannot exceed a certain value MJ at < 1 
specified by the condition <j){Md,p) = 0. 

It is instructive to compare our results with those pro- 
vided by he hydrodynamic " choking" saturation [14] and 
the conventional DVS mechanism separately. Without 
DVS (i.e. p = 0), wc reproduce the "choking" criterion 
[14] as Mj ot = 1. In the conventional hydrodynamics, 
the above condition corresponds to a subsonic gas flow 
which cannot evaluate to a supersonic one. In the oppo- 
site DVS case, (see, for example, Rcf.[15]) the current is 
saturated when the condition 1 — psM = is satisfied at 
the drain contact. We emphasize that both the mecha- 
nisms combined (as in Eq.(6)) result in a faster current 
saturation, compared with the case when each of them is 
considered separately 

The solution to Eq. (6) determines the steady-state spa- 
tial dependence of the Mach number M(rj): 
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where f{M) = M~ 4 / 3 + 2M 2 / 3 - 2pA/ s 1/3 M~ 2 / 3 . At a 
fixed friction strength 7, the Mach numbers M s and Md 
are related by 



f(M s ) - f(Md) = 



2- 
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Substituting Md = MJ at into Eq.(8), we obtain the satu- 
ration threshold 7 sa t(M s ,p), represented in Fig. 3 by dot- 
ted lines, for different values of the DVS parameter p. 
As expected, the saturation current (~ M s ) at fixed 7 
decreases as the DVS parameter p increases. 



j = WCVU = const, 
dM _ 3 7 M 7/3 
"&T~ 2 M V3 0(Af, P ) : 
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where 77 = x/l is the dimensionless coordinate, 7 = 
1/(S s t) the friction parameter, and S s the plasma wave 
velocity at the source. Then, <fi(M,p) = 1 — M 2 — psM, 
where s = S/S s = (M5/M) 1 / 3 is the dimensionless local 
plasma wave velocity, and M s is the Mach number at the 



Instability 

Wc now study the stability of the steady-state flow by 
inspecting the temporal behavior of small perturbations 
Vi(n), Ui(rj) superimposed on a steady flow with a veloc- 
ity V(rj) and channel potential U(rj). Thus, we take 

V(r),t) = V(r)) + V 1 (r ) )exp(-iu;t), 

U(r), t) = U{rj) + U!(ri)exp(-iujt). (9) 
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FIG. 1: Dimensionless frequency "^J' vs the Mach number at 
the source for the "clean limit" [1] at p; 7 = 0. The shadowed 
area corresponds to that in the instability threshold diagram 
(see Fig. 3) and shows all the frequencies that can be generated 
( lu" > ) at fixed p = 1 and < 7 < 0.32. The bold line 
A-B shows the frequency ui'{M s ) in the saturation mode (see 
the corresponding line in Fig. 3). The point B corresponds to 
MT ax {l) = 0.62. 



FIG. 2: Dimensionless increment of instability lu"1/S s vs 
Mach number at the source for the: "clean limit" [1] at 
p;7 = 0, 2DEG with DVS(p = 1) and different values of 
the friction parameter 7 = 0; 0.1; 0.2; 0.3. The bold line A- 
B shows the increment of instability in the saturation mode 
(represented by A-B line in Fig. 3). The point B corresponds 
to M™ ax {l) = 0.62. 



Equations (1),(2),(3), linearized with respect to Vi and 
Ui, yield: 

- iQ~/Vi + -^-(v vi + ui) - p(viu' + Vqu'i) + 7U1 = 0, 

— ifl'yu-i + — (uqVi + uiv ) = 0. 
ar) 

uo(0) =0,uo(l)wi(l)+«i(l)«d(l) = Q[10) 

where vi = Vi/S s ,Vq = V(r])/S s ,ui = Ui/U s ,u Q = 
U(rf)/U a . Then, lu = Cl/r = Sl-yS s /l, where ft is the 
dimensionless frequency. The complex frequency SI = 
D.' + iQ," can be determined by solving these equations. 
A positive imaginary part f2" > corresponds to insta- 
bility. 

It is instructive to introduce the current perturbation 
w = UqVi + u±vq. With the help of a simple relationship 
v = (MsM 2 ) 1 / 3 specified by Eq.(6), Eq.(10) can be re- 
written as follows: 

^(M, + „^( 2i «-l±^) + ,n, 
n + i— —— = 0, 

w'{0) = 0,w(l) = 0, (12) 

We remind that the scale of the steady state nonuni- 
formity is the channel length. The wavelengths of the 
lower-order modes of the plasma oscillations are on the 
same order of magnitude. In contrast, the higher-order 
modes have a considerably shorter scale. As stated in 
Rcf. [4] . the higher modes are more sensitive to the steady 
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FIG. 3: Diagram of instability threshold. The critical value 
of the friction parameter, 7 cr vs Mach number at source (bold 
lines) and saturation threshold 7 sa t (dotted lines) for p — 
(see Ref. [4]) and p = 1. The thin line is the asymptote 
7 cr = 2M a . The increment of instability along the line A-B is 
shown in Fig.2. The point B corresponds to M s maa: (l) = 0.62. 
The ultra-bold line represents the maximum of the instabil- 
ity threshold y m ax{M 3 ) for < p < 00. The data obtained 
by numerical simulations based on Eqs.(l,2) are shown by 
squares. Inset: FET parameters 7 vs 7* required for plasma 
wave generation. 



state nonuniformity than the fundamental and low-order 
modes. Since the increment of instability is strongly af- 
fected by the boundary-reflection coefficients, it would 
be expected that the higher-order modes are amplified 
more effectively. Based on such a qualitative argumen- 
tation, we analyze Eq.(ll) within the higher-order mode 
approximation. 

The spatial dependence of the Mach number M(rj) is 
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a monotonous function of the coordinate. Consequently, 
we can represent the current perturbation w as a func- 
tion of the Mach number w = w(M). Moreover, fol- 
lowing Ref. [4] we search for the solution in the form 
w = Fexp(g), where g = + 5/4)ln(Af). Finally, 

we use Eq.(6) and, then, change the independent variable 
from T) to M, to rewrite Eq.(ll) as follows: 



-M 4 F" 
4 
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F = 0. 



w'(M a ) = 0,w(M d ) =0. (13) 

where we now denote by F' , F" ,w' ,w" the derivatives 
with respect to M. It is noteworthy that Eq.(13) is 
exact. Solutions to this equation with given bound- 
ary conditions exist only for certain complex frequencies 
H (eigenvalues of the problem). Note the analogy be- 
tween Eq.(13) and one-dimensional Schrdingcr equation, 
in which the complex frequency f2 plays the role of energy 
in the quantum mechanics. In general, Eq.(13) can be 
solved only numerically. However, an analytical solution 
can be found using the high-order mode method [4] which 
corresponds to the quasi-classical approximation in quan- 
tum mechanics. This approach is valid when 3> 1, 
and, therefore, the imaginary part of f2 is small as com- 
pared with its real part, SI' ^ SI". We now search for a 
solution to Eq.(13) in the form F = ^= exp(— i j kdM), 
where A is a constant, and k(M) is the dimensionlcss 
wave vector. Finally, we obtain the dispersion equation 
as 



h,2 = ± 



MP 



n^/i-psM 



2V1 - psM 



(14) 



Here we keep the main corrections with respect to pow- 
ers of S7, and, therefore, neglect the term j§M 2 in 
Eq.(13). The wave vectors ki,k 2 correspond to the 
plasma waves propagating down- and upstream, respec- 
tively. Searching the solution to Eq.(13) in the form 
F = J2 ^tr ex PH JhdM), with the given bound- 



1=1.2 vl fc 'l 
ary conditions, we can easily find 



cxp U J (k 2 - k^dMJ = i£R£ 



(15) 



where R{ 



y/l-psM+M s D j 



-1 are the current am- 



plitude reflection coefficients from the source and drain 
boundaries, respectively, and R a R d is the round-trip gain 
factor. Eq.(15) allows finding of both the real and imag- 
inary parts of S7 as follows: 

O' = 3nn 1 Mfil 

M \M. 



M \M. 



We emphasize that, in the "clean limit" [1], i.e., at 
j:p = 0, Eqs.(16,17) coincide with those specified by 
Eq.(4). Using of Eq.(17) and the relationship between the 
Mach numbers at the source and drain given by Eq.(8), 
we can easily obtain the dependence of the instability 
increment u>" on Mach number at the source M s for ar- 
bitrary friction 7 and DVS strength p. Both the fric- 
tion and the DVS effect result in instability damping (see 
Fig. 2). Firstly, DVS leads to narrowing of the instability 
range. For example, at p = 1 and 7 = the instabil- 
ity exists within the range < M s < A/] nQX (l) = 0.62, 
which is narrower than that < M s < 1 in the "clean 
limit". Secondly, at certain DVS strength the enhance- 
ment of friction suppresses the instability (see the curves 
for p = 1 in Fig. 2). The region of instability narrows as 
well. Actually, in the case of p = 1 (Fig. 2), the instability 
exists when 7 < 0.32. This condition is more stringent 
as compared with the analogous 7 < 0.54 in the absence 
(p = 0) of the DVS mechanism. 

Substituting 0" = into Eq.(17) and using Eq.(8), 
we plot in Fig. 3 the instability threshold 7 cr (M s ,p) for 
two different DVS strengths p = 0; 1. In addition, the 
same figure shows by dotted lines the corresponding de- 
pendencies for the saturation threshold r y sa t(M s ,p). At 
a certain point (denoted by " A" for the p = 1 case) both 
the threshold curves merge, and, thus give the threshold 
diagram of instability, marked by the shadowed area in 
Fig. 3. The threshold diagram is confined between two 
asymptotes. One of these (thin line in Fig. 3), given by 
7c,- = 2M S , is the same for both p = and p = 1 cases. 
Actually, this is the result obtained previously in Ref. [1] 
for the case of a small friction, at which the spatial dis- 
tribution of M is only slightly nonuniform. Then, the 
threshold diagram is bounded at high M s by the satu- 
ration threshold curve 7 so t(M s ,p). Note that, within a 
part of the latter asymptote (for the p = 1 case, the bold 
line A-B in Fig. 3), the increment of instability is positive, 
as shown in Fig. 2. In what follows, we explain the above 
discrepancy by the uncertainty of our high-order mode 
method. 

Arguing that Fig. 3 represents the instability thresh- 
old, we, in addition, determine 7 cr by numerical simula- 
tions based on Eqs.(l,2) and Eq.(3), using the previously 
described numerical method [2]. The results obtained 
(shown by squares in Fig. 3) coincide with our theoretical 
curve within the error of the numerical method. As the 
numerical simulations give the threshold for the most un- 
stable modes, this coincidence with the analytically ob- 
tained results, clearly demonstrates the validity of our 
approach. 

We now present our central result, concerning the de- 
vice parameters required for a possible observation of the 
instability. First, for certain p we find the maximum of 
the instability threshold ^y m ax(p) (the corresponding de- 
pendence jmax{Ms) is represented by the ultra-bold line 
in Fig. 3). Then, the general instability condition reads 
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7 < lmax{p)- It will be recalled that p = 2*., where 
we introduce the dimensionless parameter 7* = - — - re- 
lated to the DVS effect. With the above notation, the 
instability criterion yields 

7*>7™L(7)7, (18) 

where n j^ lx {x) is the inverse function. In the inset in 
Fig. 3, we plot the diagram specified by Eq.(18). The val- 
ues of the FET parameters, confined by the bold curve, 
correspond to an instability. Another condition is related 
to the applicability of the hydrodynamic model suggested 
in Ref. [1] . In actual fact, 2D electrons can be considered a 
fluid when the mean free path A ~ y/N, associated with 
electron-electron collisions, is shorter than the channel 
length. This requirement can be easily transformed to 

The final question remaining to be resolved is whether 
the finite scattering and the DVS effect can influence the 
proper and the higher-order-mode frequencies of the in- 
stability. With the help of Eq.(16), we calculate the nth- 
order- mode frequency uj' n (M s ,p, 7). At a fixed p = 1 
in Fig. I, we plot the shadowed area (which corresponds 
to that in Fig. 3) which shows all the possible plasma 
frequencies when the instability increment is positive 
51" > 0. At a fixed current j ~ M s and a finite p, the 
proper frequency and its harmonics may be much lower 
than those expected in "clean limit" (see Eq.(4)). The 
difference between two cases becomes more pronounced 
in the saturation mode, at M 8 < M™ ax . We conclude 
that the DVS effect results in a strong depression of the 
plasma wave generation frequency. 

We are now going to justify the applicability of the 
high-ordcr-mode approach[4], developed in the present 
paper. We first compare the data presented in Figs. 1, 2, 
and, then conclude that our basic assumption f2' > 51" 
is, surprisingly, valid for the proper plasma wave mode 
n = 1 as well. This is, however, not the case in the 
saturation mode, at f2 — ► 0, which seems to be the reason 
for a minor deviation between the theory and the results 
of our numerical simulations (see threshold diagram in 
Fig.3). 

Let us represent the instability threshold shown in 
Fig.3 in current-voltage ( I-V ) plane of the transistor. 
In terms of the Mach number formalism, the current and 
the source-drain voltage swing yield 

— = 7- 3 M„ ^ = 7 - 2 (l - {M 8 /M d f/% (19) 
Jo U 

where we use the units j = WCm(l/T) 3 /e and Uo = 
m(l/T) 2 /e, introduced in Ref. [14]. For a fixed M s ,7,p, 
we can calculate Md specified by Eq.(8) and then 
use Eq.(19) to find the steady-state I-V characteristics 
j(U s d)- In Fig. 4 we plot I-V characteristic for certain 
values p = 0; 1. In addition, the instability threshold 
curves j cr (M s ,p) can be transformed to those j C r{U s d) 
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FIG. 4: Instability region in the I-V characteristic. The I-V 
curves are plotted for 7 = 0.3, where the upper(lower) curve 
corresponds to p = and p = 1, respectively. For p = and 
p = 1 cases the bold (dotted) line (corresponding to those in 
Fig.3) indicates the instability (saturation) threshold respec- 
tively. The unstable parts of the current-voltage characteris- 
tics are indicated by thin lines. The line A-B corresponds to 
that shown in Fig.3 

in the I-V plane. Furthermore, the dotted curves (cor- 
responding to 'y sat (M s ,p) in Fig.3) relate the saturation 
current to saturation source-drain voltage. For the p = 1 
case, both the threshold curves, which merge at point A, 
confine the instability region in the I-V plane. The parts 
of the current-voltage characteristics lying within this re- 
gion correspond to unstable states (thin lines in Fig. 4). 
For 7 > 0.32(p = 1) and 7 > 0.54(p = 0), the current- 
voltage characteristics are always stable. Note that the 
saturation current at finite p is lower than that that in 
absence of the DVS effect. 



Theory vs experiment 

We now provide a detailed analysis of the data 
obtained in a THz-cmission experiment[12] for an 
ultrashort-channel InGaAs-bascd FET. In this device, 
the source-drain separation of 1.3^m was made much 
longer than the gate length I = 60nm. Thus, the total 
series resistance R associated with the gate-source and 
gate-drain access regions may affect the I-V characteris- 
tics of the device. Indeed, the voltage bias between the 
source and drain ends of the gated part of the channel 
U s d is related to the external (hereafter marked by prime) 
source-to-drain terminal voltage by U' sd = U s d + jR- 
Then, for the actual case of a symmetric FET (see inset 
in Fig. 5) the swing voltage U s is related to the exter- 
nal gate-to-source voltage by U s = U' gs — Ut — jR/2. 
We further use the simplest channel pinch-off model, 
for which the I-V dependence of the device is given by 
j = l3U sd {U s - U sd 1 '2), where = WCfi/l is the I- 
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TABLE I: Experiment vs theory 



U' gs (mV) 



FIG. 5: I-V transfer characteristic ( bold line ) at U' ad = 5mV, 
reproduced from Ref. [12]. Dotted line represents the fit with 
Eq.(20). Inset: schematic of the short-channel FET. 



V steepness. Assuming that the actually measured I-V 
characteristic is expressed in terms of terminal voltages 



U' d , we derive I-V dependence as follows 



3=P 



1 + 0R(UI-U'J2Y 



(20) 



where U' s = U' gs — Ut- We now analyze the I-V transfer 
characteristics measured [12] at a fixed source-drain bias 
U' sd = 5mV. According to Eq.(20), if the applied gate-to- 
source terminal voltage U' gg exceeds a certain threshold 
value Ut + U' sd /2, the current is nonzero. Then, at high 
U' gs , the current saturates, being limited only by the only 
series resistance, and, hence j = U' sd /R. In Fig. 5, we 
reproduce the observed FET transfer characteristic (see 
inset in Fig.l of [12]) and then put on the same plot our 
best-fit data. Note that Eq.(20) is, in fact, valid only 
in the gradual Shockley approximation, and, therefore, 
the subthreshold mode U' — Ut < is not analyzed in 
Fig. 5. The device parameters obtained are as follows: 
(3 = 0M^y,R = 140, and U T = -147mV. For e = 12.7, 
gate-to-channel separation d = 17nm, and gate width 
W = 50/mi, the FET mobility is fx = 1200cm 2 /Vs. The 
corresponding momentum relaxation time is r = 2.9 x 
l(T 14 s. 

The next step of our analysis is concerned the FET I- 
V characteristic [12] at a fixed gate-to-source voltage (in 
the experiment U' gs = 0, and, therefore, U' s = —Ut)- At 
a fixed U' s the current is linear in the source-to-drain bias 

j = R U ^ R , where R c h = {pu's)" 1 denotes the channel 
resistance. At a certain source-to-drain bias, the device 
current saturates (channel pinch-off) when dj/dU' sd = 0. 
Both the saturation current and voltage are given by 
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where r = R/R c h- Without the series resistance R, 
Eq.(21) reproduces the conventional channel pinch-off 
saturation at j sa t = PU^/2, Uffi = U s . By contrast, 
at a nonzero R, the saturation current decreases and the 
saturation voltage increases as compared with the R = 
case. For the I-V characteristic reported in [12], we have 
U' s = 147mV, and, therefore, the channel resistance is 
R ch = 110. The total resistance R + R c h = 14+ 11 = 25 
is comparable with the value of 290, obtained in [12] 
from the linear part of the I-V curve. Then, we esti- 
mate both the saturation current j sat = 4.4mA and volt- 
age U' sd \ sat = 180mV, which are, however, somewhat low 
than those obtained in the experiment. 

We now discuss the most intriguing result [12] concern- 
ing the irradiation / ~ 0.4THz, observed at certain cur- 
rent j = 4.5mA and voltage U' sd = 200mV. The above 
emission threshold corresponds to U sd = U' sd — jR = 
70mV and effective voltage swing U s = U' gs — UT—jR/2 = 
116mV. In contrast to Ref. [12], we conclude that the ob- 
served emission threshold precedes the saturation mode 
of the transistor because U s d < U s . Our conclusion can 
be further confirmed by analyzing the actual shape of 
the I-V characteristic reported in Ref. [12]. In general, 
the output drain conductance vanishes in the saturation 
mode, i.e. d ^, = 0, irrespective of the actual mechanism 
(channel pinch-off, DVS mechanism based on smooth 
t(E) dependence, "chocking"). Obviously, this is not 
the case for the THz-emission threshold point reported 
in Ref. [12]. 

We now verify whether the device parameters 7, 7* 
match the limits( see Fig. 4, inset ) required for plasma 
wave instability to occur. For m = 0.042mo, we find the 
plasma wave velocity S s = ^JeU s /m = 7.0 x 10 5 m/s 



and the 2D density N s 



e n eU B 
de 



3.5 x 10 1 



We estimate the friction strength as 7 = -J^- = 3.0. 
Then, with the textbook value for the saturation ve- 
locity V sa t = 3.2 x 10 5 m/s, we obtain p = 2.2, hence 
7* = 6.6. Both the parameters 7,7* exceed the limits 
( see Fig. 4, inset ) required for instability. Hence, the 
observed sub-THz emission [12] cannot be attributed to 
the plasma wave instability scenario. Moreover, our step- 
by-step analysis (see Table 1) of the data for a 1.5/xm- 
channel GaN FET [12] shows that the GHz-cmission ob- 



7 



served does not originate from the " shallow- water" mech- 
anism, cither. 



Conclusion 

We have shown that both the scattering and the DVS 
effect depress the frequency and the instability incre- 
ment. The experimental results reported in Refs.[12],[12] 
for sub-THz emission in ultra-short channel FETs can- 
not be explained in terms of the plasma instability sce- 
nario [1]. 
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